Rate of steady— state reconnection in an incompressible plasma 



Nikolai V. Erkaev 

Institute of Computational Modelling, Russian Academy of Sciences, 660036 Krasnoyarsk 36, Russia 

Vladimir S. Semenov, Ilya V. Alexeev 
Institute of Physics, University of St. Petersburg, St. Petergof, 198504 • Russia 

Helfried K. Biernat 

Space Research Institute, Austrian Academy of Sciences, Schmiedlstrasse 6, A-8042 Graz, Austria 

The reconnection rate is obtained for the simplest case of 2D symmetric reconnection in an 
incompressible plasma. In the short note (Erkaev et al.El), the reconnection rate is found by matching 
the outer Petschek solution and the inner diffusion region solution. Here the details of the numerical 
simulation of the diffusion region are presented and the asymptotic procedure which is used for 
deriving the reconnection rate is described. The reconnection rate is obtained as a decreasing 
function of the diffusion region length. For a sufficiently large diffusion region scale, the reconnection 
rate becomes close to that obtained in the Sweet-Parker solution with the inverse square root 
dependence on the magnetic Reynolds number Rcm, determined for the global size of the current 
sheet. On the other hand, for a small diffusion region length scale, the reconnection rate turns out 
to be very similar to that obtained in the Petschek model with a logarithmic dependence on the 
magnetic Reynolds number Rcm- This means that the Petschek regime seems to be possible only in 
the case of a strongly localized conductivity corresponding to a small scale of the diffusion region. 

I. INTRODUCTION 

Magnetic reconnection is a physical process in plasmas which changes a magnetic field topology and releases stored 
magnetic energv. It is one of the central concerns in astrophysical, solar, space, fusion and laboratory plasmas (e.g., 
Honeaj; PriestEl). 

A key question arising in the reconnection theory is that of the reconnection rate. So far there are two different 
magnetohydrodynamic (MHD) models of reconnection based on the Sweet-Parker (pure diffusion) (see ParkeiO; 
SweetO) and the Petschek (slow shock energy conversion) (see PetschekH) approaches. These models propose two 
different estimations of the reconnection rate e: The Sweet-Parker model predicts e ~ l/y/Rem, and the Petschek 
model gives e ~ l/lni?em, where 

Rdni = 5 (1) 

is the global magnetic Reynolds number based on the half-length of a current layer L, the Alfven velocity Va, and 
the resistivity of the plasma -q. For cosmic plasmas, magnetic Reynolds numbers usually are very large, therefore the 
Petschek regime seems to be much more effective. However, since the Petschek reconnection model was proposed, it 
is not clear what conditions are necessary to realiza this regime. 

It is a fact that numerical simulations (BiskampQ; Scholem) carried out for a constant resistivity were not able to 
reproduce the solution of Petschek type, instead, they were rather in favour of the Sweet-Parker solution. Laboratory 
experiments also seem to observe the Sweet-Parker regime of reconnection (Ji et al.u). 

On the othei^iand, if nonuniform resistivity is localized to a small region, the results of numerical simulations 
(ScholerS; Ugail3) clearly show Petschek-type reconnection with pronounced slow shocks. For the Petschek regime, 
there are two physically different regions: A small diffusion region, where dissipation is important, is surrounded by a 
large convective zone where the plasma can be considered as ideal and dissipationless. The problem is very complicated 
and thus it does not seem realistic to obtain an analytical solution which is valid for both regions simultaneously. 
To simplify this problem, we seek solutions separately, in the diffusion region and in the convective zone. For the 
later, a solution can be obtained analytically as an asymptotic series with respect to a small reconnection rate. For 
the diffusion region, it is impossible to find an analytical solution, and hence it has to be obtained numerically. In 
this semi-analytical approach, we have to combine the numerical solution for the diffusion region and the Petschek 
analytical solution for the convective region. The latter can be done by different methods, which lead to absolutely 
identical results for the reconnection rate estimation. The estimation obtained by Erkaev et alJ^I is based on asymptotic 
matching of the diffusion region and convective zone solutions. In our present work, we use another way based on a 
regularized convective region solution, which seems to be rather clear and very close to the original Petschek method. 
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In this paper we give a detailed description of the numerical solution for the diffusion region, and derive the estimation 
for the reconnection rate. 

This paper is organized as follows: In Sections II and III, we start with the steady-state MHD equations and present 
the Petschek solution. The diffusion region scaling and boundary layer equations arc introduced in Section IV. The 
numerical algorithm and the results of the calculations are described in Sections V and VI. The reconnection rate 
is derived in Section VII, whereas Section VIII is devoted to the summary and discussion. Mathematical details are 
described in the Appendix. 



II. MHD EQUATIONS 

In the problem under consideration, the plasma is governed by the resistive steady-state MHD system of equations 

p(vV)v = -VP + -^(B-V)B, (2) 

1 c 

E+ -(v X B) = — r/(a;,y)curlB, (3) 
c 47r 

V • B = 0, V • V 0, (4) 

where p is a mass density, P is the total pressure, P ^ p + B'^/Stt, and Rem is the global magnetic Reynolds number 
based on the maximal value of the resistivity rjmax ■ 

Outside of the diffusion region, in the so-called convection zone, dissipation is not important any longer, and we 
can use the ideal system of MHD equations in the limit Rem — *■ oo. 

In an incompressible plasma the following relations have to be satisfied at the shock front 



{Bn} = 0, (5) 

{Vn} = 0, (6) 

{P} = 0, (7) 

^BnBt - pvn^rtj = 0, (8) 

{B„vt - w„Bt} = 0, (9) 



where the subscripts n and t denote components normal and tangential to the shock front. 



III. PETSCHEK SOLUTION 

The Petschek solution, which is valid in the convection region, can be presented as follows (Petscheki, for details see 
Vasyliunaaiil) . We use coordinates x, y, which are directed along the current sheet and in the perpendicular direction, 
respectively. The solution is completely determined by the following parameters: Quantity L which is the halflength 
of the current sheet, Vq is the plasma inflow velocity, and Bq is the initial magnetic field. The solution is presented in 
the form of asymptotic series with respect to the small parameter which is known as the reconnection rate 

Here Eq is the electric field which is constant in the 2D case under consideration, and Ea — ^VaBq is the Alfven 
electric field. 



Inflow region: 



Vj; = 0, Vy^-eVA, (11) 



— , , By — arctan — 

+2/2 TT y 



Bx = Bq — In — ^ ^ , By — — - — arctan — . (12) 



Outflow region: 



Vx = Va, Vy = 0, (13) 
B, = 0, By = eBo. (14) 
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The equation for the shock in the first quadrant is 

y^ex. (15) 
It can be shown that sHghtly outside of the shock from the inflow side 

5.(^,0) = I _2eBo:, < ^^^^ 

Expressions (|Tl[-p^ are asymptotic solutions with respect to e (zero and first order terms in tlic inflow region and 
only zero order term in the outflow region) of the ideal MHD system of equations (^-|^) and the Rankine-Hugoniot 
shock relations (|[^). 

Petschek did not obtain a solution in the diffusion region, instead, he estimated the maximum reconnection rate as 
l/ln_Rem using some simple physical suggestions. Generally speaking, this implies that the Petschek model gives any 
reconnection rate from the Sweet-Parker value 1 / V Rem up to l/lni?em, and for a long time, it was unclear whether 
Petschek reconnection faster than Sweet-Parker reconnection is possible. This problem can be solved by combining 
the analytical Petschek solution ( |lT[jl^ ) and the numerical model of the diffusion region. 



IV. DIFFUSION REGION SCALING 



The next step is to find a numerical solution for the diffusion region. But first we have to obtain the boundary 
layer MHD equations suitable for the diffusion region. 

To this end we renormalize the MHD equations to new scales Ed, VacI, EacI — BdVAd/c, Pd, where all quantities 
are supposed to be taken at the upper boundary of the diffusion region, 

x'^x/l^, y' = y/l^, B' = B/Bd, V = w/VAd, P' = P/Pd, (17) 

where Z,, is the characteristic length of the resistivity variation. The diffusion region length scale Id (see Figure 1) 
obtained from our numerical results (Section VII) is of order of the scale 

The convective electric field — v x B/c is zero in the centre of the diffusion region x = y = Q where v = B = 0, 
and then increases to the constant value Eq at the boundary of the convection zone. This type of behaviour of the 
convective electric field is reasonable to be used for the definition of the size of the diffusion region which is one of 
the most important parameters of the problem. Namely, the length scale of the diffusion region is determined as the 
distance between the origin x = 0, y = and the boundary where the convective electric field reaches its asymptotic 
value i?o, or better to say, some level, for example, Q.QEq. 

In the diffusion region where dissipation is essential, we adopt the dissipative MHD equations with the magnetic 
Reynolds number 

Re'd - > (18) 

C f]max 

and the normalized electric field E' = Ec/ {VAdBd) — e', where e' is a local reconnection rate at the diffusion region 
boundary. These electric field and local reconnection rate are not known. They are to be obtained from the numerical 
solution for the diffusion region. I— . 

The scaling for the diffusion region is similar to that for the Prandtl viscous layer (Landau and LifschitzllHl) and 
corresponds exactly to the Sweet-Parker one 

/ TDi ..I ni _ /n/i \ 

(19) 




Consequently, the new boundary layer variables are as follows 

x^x', Bx^B'^, Vcc^v'^, P = P', 

y^y'^/lK. By^B'y^/R^^, Vy^v'y^/EZ^, e = e'/R^. ^ > 

The diffusion region Reynolds number Re'^ is certainly smaller then the global Reynolds number Rem, but still it 
is supposed to be Re'^ S> 1. Therefore, in zero-order with respect to l/i?e'^, the boundary layer equations turn out 
to be 
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dvx , dvx , dvx ^ dBx ^ dB^ dP 
dt ^ dx ^ dy ^ dx ^ dy dx 
dP 

"^i = 0, (25) 



dx dy 

= 0, (26) 



dvx dvy 



dx dy 

where f}{x^y) is the normahzed resistivity of the plasma with the maximum value to be 1, ^ is a small parameter, 
/I = l/i?e^. The small terms which include /i at the right sides of the induction equations are necessary for numerical 
stability of the calculations. 

It can be seen from equation (|2^) that the total pressure is constant across the diffusion region. This is a general 
feature of a boundary layer approximation. Hence, the total pressure is defined inside the diffusion region by values 



at the boundary, and for the boundary layer equations (|2l|-|26|), the total pressure can be considered to be a given 
function of x, e.g., P{x). 

As it was pointed out, the appropriate exact solutions of the boundary layer equations (pl|-p6|) are unknown even in 
the steady-state case, therefore, we have to solve the problem numerically. Although we have to obtain a steady-state 
solution, from the point of view of simulation, it is advantageous to use a relaxation method and to solve numerically 
the unsteady system of the boundary layer MHD equation s (2l|-2^ ). 



It is important to note that in the subset of equations (p3~|25), only two equations are independent. In principle, 
we can determine the normal component from the induction equation ( p^ ) or from the equation ( p5| ) providing the 
magnetic flux conservation. From the mathematical point of view, they are equivalent. In our numerical solution, 
we use equation ( p5| ) to determine the By component in the internal grid points, and the equation ( |2^ ) is used as a 
boundary condition at the lower boundary. 



V. NUMERICAL ALGORITHM 



Starting with an initial MHD configuration under fixed boundary conditions, we look for the convergence of the 
time-dependent solution to a steady state. To avoid additional numerical diffusion, we do not use a flux function and 
a magnetic potential. The normalized total pressure is chosen to be 1. ^ |__. 

The distribution of the resistivity 77 = rjmaxVix, y) is traditional (Scholeia; UgaillEl) 

7y(i,y) = de(-^-^'-^='^')+/, (27) 

with d + / = 1. Setting d = .95 and / = .05 we can model a case of localized resistivity, for c? = and / = 1 the 
resistivity is uniform. 

As the initial configuration, we choose a current sheet with a linear profile of the magnetic field Bx — y, By — 0. 
The velocity components are assumed to be equal to zero at the initial moment, Vx = 0, Vy = 0. 

To solve the MHD system numerically, we use a two step conservative finite difference numerical scheme with a 
rectangular grid 145 x 100 in the first quadrant. From a time level (n), we calculate the parameters on the next time 
level (n + 1) in two steps. In the first step (n + 1/2), diffusion is switched off, and we calculate the parameter at the 
intermediate points (n + 1/2) using the equations in characteristic form. This is similar to the approach used in the 
Godunov method. In the second step, we calculate the parameters at the next time level (n + 1) using the equations 
in conservative form and taking into account the diffusion terms approximated in implicit form. 

The details of the numerical algorithm are the following. The Bx component is found from the x-component of the 
induction equation 



d .^dBx\ 



4 



d / ..dBy 

TTz ri{x,y) 



ay \ OX 



(28) 



where the parameters 

= {B,j:Vy - "i4i?y)"t+i/2 (29) 

are determined by the method of characteristics on the level n+ 1/2. This implies that at the beginning {n n+1/2) 
diffusion is switched off, and only convection acts, and then for given convection, diffusion is switched on, and is 
calculated on the level n. The normal magnetic field component By is determined from the equation V • B = 
The velocity component Vx is found from the x-component of the momentum equation ( pT| ) , 

{iy-)lt' - iy-)7,k]/r + {Qy^,k+l/2 - QyrM-l/2r^'^Vhy + (Q„+i/2,fc - Q y^-l/2.kr^'^V hx = 0, (30) 

where 

Qytk+l/2 ^ i'^xVy - BxBy).J^^{^^ , (31) 

Qn+1/2 _ ,y2 _ n2\n+l/2 /oo) 

Here, the parameters ()"^^(/9 are determined by the method of characteristics on the level n + 1/2 simultaneously 

with the calculation of Bx- The velocity component Vy is determined from the equation divV = 0. 

The boundary conditions are as follows: 
At the upper (inflow) boundary, the tangential magnetic field component is assumed to be constant, B^ =1 and the 
tangential velocity component vanishes Vx — Q- 

At the left boundary we have the symmetry conditions, dBx/dx = 0, By ~ 0, Vx = 0. 

At the right boundary we hold free conditions suitable for a uniform flow in the outflow region, dBy/dx — 0, 
dVy/dx = 0. 

At the lower boundary (y = 0) there is the symmetry condition for the tangential magnetic field component, Bx — 0, 
and the non-flow condition for the normal velocity component, Vy = 0. At this boundary, the normal component of 
the magnetic field By is obtained from the induction equation (^ij) on the line y — 0, 

dBy d , d ( ,^ .,dBx\ d ( ,^ ~.dBy\ 

-i- + = -3=x j j ■ 

The small parameter /i ^ .1 — .2 is used here to regularize the numerical scheme for the unsteady system of the 
boundary layer MHD equations (pT|-p^), which is an ill-posed problem in our case. 

The size of the computational domain is chosen to be much larger than the diffusion region size Id-, and also much 
less than the global size L. At the inflow boundary we do not fix the normal components of the magnetic field and 
velocity, and thus we do not impose a reconnection rate and an electric field in the diffusion region from the very 
beginning. The latter has to be found from the numerical solution self-consistently. 



VI. RESULTS OF THE NUMERICAL SIMULATION 



To estimate the convergence of the time-dependent solution to a steady state for each n-th time step, we use 
the following criteria, maxdV!^" — Vx~^\) I {l^t\Vxmax\) < 10^®- In the 2D steady state the total (convective plus 
dissipative) electric field must be constant, and it is so in our simulations (see Figures 2, and 3) besides of small 
perturbations near the outflow boundary due to some reflections, although we apply free boundary conditions. 

Let us discuss the result of our simulations. For the case of localized resistivity, the system reaches the Petschek 
steady state with clear asymptotic behaviour (see Figure 2): V^j; ^ 1 in the outflow region; V^j, ^ e at the inflow 
boundary; B^ decreases from 1 to at the shock transition; By ^ e in the outflow region; and By — > 2e from the 
inflow side of the shock (compare with the Petschek solution (pT| p^). 

There is a well pronounced slow shock, as can be seen in the behaviour of all MHD parameters, but in particular 
in the distribution of the current density. The normalized electric field (reconnection rate) turns out to be £ ~ 0.7. It 
is important to note that the numerical results do not depend on the size of calculation box. 



5 



On the other hand, for the case of homogeneous resistivity, the system reaches the Sweet-Parker state (see Figure 
3), even if the Petschek solution is used as initial configuration (see also Scholem; UgaillJ; Uzdensky and Kulsruc£3). 
This seems to imply that Petschek-type reconnection is possible only if the resistivity of the plasma is localized to a 
small region, whereas for constant resistivity, the Sweet-Parker regime is realized (Erkaev et al.El). 

The size of the diffusion region layer Id is defined as its length along the x axis where the convective electric field at 
the lower boundary (y — 0) Ec — —v^By is less in absolute value than some level of the total electric field (say 0.9e). 
For the case of a localized resistivity, Id practically coincides with the scale of the inhomogeneity of the resistivity 
Iri when the maximum of resistivity is much larger then the background resistivity. Therefore hereafter we consider 

For the case of uniform resistivity, the plasma is accelerated very slowly, and there is no obvious definition for the 
scale length of the diffusion region. Diffusion is important everywhere for the pure Sweet-Parker regime, and for the 
Petschek asymptotic solution there is left no room. Therefore, the solution does not converge to the Petschek solution, 
not only at the right hand boundary but everywhere. In this case, the solution will depend on the calculation box size 
because it does not have any other scale. Hence, the constant resistivity solution can not be matched to the Petschek 
solution. 

Nevertheless, the Sweet-Parker regime is still important also for the Petschek solution, because in the nearest 
vicinity of the reconnection line, where the resistivity can be considered to be constant, the diffusion region structure 
is similar to the Sweet-Parker case. Besides, and this is even more important, the scaling for the diffusion region is 
exactly the Sweet-Parker one (^, or, better to say, the Prandtl scahng. 



VII. RECONNECTION RATE 



To find a relationship between the reconnection rate and dissipation we need first of all an estimation of magnetic 
field at the boundary of the diffusion region Ed- To this end we can not use the Petschek solution (|l2[) because the E^ 
component diverges at the origin Ex — > — oo, when r — \/x^\-y^ — > 0. This singularity is a consequence of the fact 
that dissipation actually has not been taken into account for the Petschek solution. Formally it follows from the jump 
at the origin of the Ey component of the magnetic field ( p^ . Dissipation evidently leads to smooth behaviour of the 
magnetic field in the diffusion region, and then no singularities are possible. To illustrate this we consider a model 
distribution of the Ep.{x,0) component with linearly smoothed boundary condition at the interval {—ld,ld) similar to 
the original PetscheUS consideration 

r ±2eEQ L>\x\ > Id 
E^ix^O)^! 2eEof^ \x\ < h (34) 
[ |a;| > L . 

The Ex(x^ y) component of the magnetic field in the inflow region can be found from the Poisson integral, 

1 y EP{x',Q){x'-: 

IT J (x' — xY + 

— oo 

n 2eBo / X {x- Idf + y2 , ^ x-ld ^ x + U 

Eq ; — [ 2ld+ - m —5 + y arctan y arctan 

-Kid V 2 [x + ldY + y^ y y 



1 f EP(x',0)(x' -x) , 
Exix,y) = E„-- / , 2 dx' 



sBof {y^ + {L^xf){v^ + {L + xf) 
T {y^ + {ld-xf){y^ + {ld + xf) 



This solution does not have a singularity at the origin any more, and tends to the Petschek solution (|T^ outside the 
diffusion region. We can simplify equation (pq) at the origin 



5,(0,0) -Bo ^In- ^. (36) 

The first term on the right hand side of this equation is of the oder of 0(1), the third one is of 0(e), but the second 
term consists of a large parameter In times the small parameter e. Thus we assume the following relations between 
the parameters 
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1 > eln — > e. 



(37) 



So far we considered only a model distribution of the By{x,0) (|3^ ) along the current sheet but it turns out that 
Bx{0,0) does not depend on the actual distribution of the By component inside the diffusion region up to 0{e). This 
implies that we can extend equation (^) to the general case. 

Let us consider the Poisson integral with the actual distribution of the By(x, 0) component using the model boundary 
condition By{x,0) (34) for regularization 



5,(0,0) = Bo-- 



+ 00 

1 fByix',01^^, 



+ 00 



Bn-- 



1 f (By ix',0) ~ BPjx', 0) + BPjx', 0)) 



+ 00 



^ 4eBo^ L 4eBo , 1 f {By{x' ,0) - B^^x' ,0)) 
Bo in \ — / ; — ax 

TT U IT TT 



Bo~^ln^-CeBo, 

TT Id 



(38) 



where C = const includes both, the contribution from and the contribution from the non-singular integral in 
the third line of this equation. The main difficulty for the estimation of this integral is that near the diffusion region, 
the local Pctschck solution reproduced in our simulation, seems to be different from the global one because e' > e 
and Bd < Bq. The local Petschek solution has asymptotically By{x/ld) — > 2e'Bd when x/ld — * oo which seems to 
be different from the condition By {x/ld) — > "^eBq used in (|3^). However, as it is shown in Appendix, the difference 
O(e') — 0(e) is of the order of e rather then 0(eln |^) (see Appendix). This allows us to estimate the integral ( |3^ ) as 
a quantity of order e which is much smaller than the main term ^ e In ^ . 

The diffusion region is small Id L and for the boundary condition for the diffusion region Bd we can use the 
magnetic field at the origin i?a;(0,0). Using the relation (^Sf), we find the magnetic field strength at the diffusion 
region boundary 



4e L 
Bd = Bo{l~-\n-). 

TT id 



(39) 



Now everything is ready to determine the reconnection rate. The electric field must be constant in the whole inflow 
region, hence 



VdBd ^ 



Wo -Bo, 
- ^Bl, 



(40) 
(41) 



where the definition of the reconnection rates e' 
(see scaling (EOh) we obtain 



Vd/Bd, £ = vq/Bq are used. Bearing in mind that e' ~ e/^/Re'^ 



eBd^l^^eB^r^^ 

y C T]max 



(42) 



Substituting Bd from equation (p9[), we determine finally the following equation for the reconnection rate e 

£(l_l£ln^)3/2=e/^, 

TT Id 



(43) 



where the magnetic Reynolds number Red — 4:TTVAld/ {c^ilmax) is based on the global Alfven velocity and the half 
length of the diffusion region Id- The internal reconnection rate e has to be found from the simulation of the diffusion 
region problem. 

For small e In ^ there is an analytical expression 
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(44) 



Here e is an internal reconnection rate, determined from the numerical solution, which is e ~ 0.7 for the Petschek 
type solution. 

In the Appendix it is also shown that the global Petschek solution with second order corrections tends to the 
asymptotic of the diffusion region solution for x Id- 

It is interesting that for the derivation of the final result (^, the only value which has been actually used is 
the internal reconnection rate e obtained from the numerical solution, and the asymptotic behaviour (^4|) . The actual 
distribution of the By component along the upper boundary of the diffusion region does not contribute at all (besides 
of the asymptotic behaviour (|3^)) in zero-order approximation considered above. Of course, from the mathematical 
point of view, it is important that the diffusion region solution exists and has the Petschek-like asymptotic behaviour 
(pl]-p^) . Therefore, the asymptotic behaviour ( |34| ) plays the key role in the derivation of the reconnection rate and 
this question needs to be clarified in more detail. 



VIII. DISCUSSION 

Equations (|4^ , ^4|) give the unique reconnection rate for known parameters of the current sheet L, Bq, Va, ?7, Id- 
Let us fix now the lengths L and start to vary Id assuming ~ Id- It is clear that for small Id, the Petschek term 
becomes large, whereas for big Id, the Sweet-Parker term is dominant. The behavior of the implicit function e{ld/ L) 
given by ( |4^ ) is non monotonic. There exists a length Id corresponding to a maximum value of the reconnection rate. 
This maximal reconnection rate is a function of the magnetic Reynolds number given in an implicit form 



4(A + ln(i?e„,/e)) ' ' ' 

where A is the constant A = 3 — 21n(e') — 31n(12/7r) = —0.31. Here Rem is the Reynolds number determined for 
the global scale and the maximal resistivity i?e,„ — '^ttVaL / {c^rjjnax) - This result can be interpreted as follows. In 
the case of a large global Reynolds number, for fixed values of the maximum resistivity and the global scale L, the 
reconnection rate and the corresponding intensity of energy conversion reach their maxima when the diffusion region 
length scale and also the conductivity length scale are much smaller than L- This maximum value of the reconnection 
rate is a logarithmic function of the global Reynolds number which is similar to that estimated by Petschek. This fact 
contradicts to the usual electrotechnical intuition. For example, to get maximum heating from a rheostat (resistor), 
we need to switch on the whole length, to increase Id, as oppose to the progress of reconnection. It is a fact that 
the energy release in the course of the reconnection process takes place not only in the form of Joule heating in the 
diffusion region and at the shock fronts, but also in the form of plasma acceleration. 

By increasing the conductivity length scale and the corresponding diffusion region length scale, the reconnection 
rate decreases substantially, becoming more close to that of the Sweet-Parker regime. 

We have to emphasize once more that the case of constant resistivity is not described by equation (^), because 
there is no clear scale of the diffusion region, no clear Petschek-type asymptotic behaviour, and therefore it can not 
be matched with the Petschek solution. 

The appearance of strongly localized resistivity is often the relevant case in space plasma applications, but for 
laboratory experiments, where the size of a device is relatively small, the Petschek regime can hardly be expected. 

One of the main difRculties of the diffusive-like theories of reconnection such as the Sweet-Parker mechanism 
(Sweet, 1958, Parker, 1963), and the tearing instability (Galeev et al., 1986) is that the efficiency of the process turns 
out to be of the order of Re~" where usually < a < 1. For example, for the Sweet-Parker regime, a — 1/2. 
In cosmic plasmas the magnetic Reynolds number is often very large because of the large scale, high velocity and 
high conductivity. Hence, the efficiency of pure dissipative processes is rather poor. The Petschek mechanism of fast 
reconnection is much more effective due to the logarithmic dependence of the reconnection rate on scale (42). In the 
Petschek model, MHD waves play the dominant role and the logarithmic dependence is the contribution of the waves 
to the efficiency of the process. 

In this paper, we studied reconnection for a strongly localized resistivity with a large ratio of the maximal and 
background resistivity ( 20). A crucial parameter for the reconnection rate is the diffusion region length which is 
obtained to be approximately equal to the length scale of the resistivity. An interesting question for future study is 
the dependence of the diffusion region length as well as the electric field on the amplitude of the resistivity variation. 
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IX. APPENDIX 



So, we have to clarify the problem concerning the asymptotic behaviour By{x/ld) — ^ 'i.s'B^ when x/ld — *■ oo, 
estimate the integral, and to prove that the global Petschek solution tends to the local one if we take into account 
all necessary terms. Originally, Petschek (1964) considered the reconnection problem using as a small parameter the 
reconnection rate e. He obtained the solution ([ll[-p^, taking into account only zero and first order terms in the inflow 
region, and zero order terms in the outflow region. But there is the possibility to extend this solution with higher 
order terms (Pudovkin and Semenov, 1985). In order to do this we have to present each component of the MHD state 
vector U (inflow region), U (outflow region), S (shock front ) as an asymptotic series with respect to the reconnection 
rate e 

C/-C/(")+eC/(i)+£2c/(2) + ... (46) 
;7 = {7(")+eC/(i)+e2c/(2) + ... (47) 
5 = + (48) 

The terms of the series (^^psf) can be obtained step by step using the MHD equations (|[^) and the shock boundary 
conditions according to the following scheme, 

t/(o) 4 ,7(0) 4 5(0) 4 t^d) 4 ^(1) 4 ^(1) 4 ^ ^ ^ (49) 

Here U'^'^^ is the initial vector, and each next term is determined via solving the reduced MHD system with boundary 
condition provided by the previous step. 

For example, the original Petschek solution (|l^^6|) corresponds to the first three steps of this scheme. The first 
step is trivial, because no shock front is yet possible. In the next step, the outflow region solution of zero order allows 
to impose a boundary condition problem for the inflow region solution in first order, and so on. 

Proceeding according to this scheme up to the step 5, we obtain the following extended Petschek solution. 

Inflow region: 

5, - So - -So In J-^ (50) 

By — — Bq arctan — (51) 
TT y 

As X 

Vx — — Vq arctan — (52) 
TT y 

Vy = -Vo - -Vo In Jl-^ (53) 



Outflow region: 



Bx^-Bo\n^, (54) 
TT X — y 

By=eB„ So In , / , (55) 

Vx^Va^ In—— , (56) 



TT ALx 

x + y ^ 

X — y X 



where y = yje. 

Shock front equation: 



/ X \ 

— 2xln-+x . (58) 

TT V L / 



Finally it is possible to find the y-component of the magnetic field By{x) at the shock which has been used in 
deriving the reconnection rate up to second order. 
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= 2Boe(l-^(ln| + 3)). (59) 



X ^ Id- From equations (39, ttlj) it follows that 



Using the extended Petschek solution (|5C-59), we can prove now that the global solution tends to the local one at 



e' = eBl(l + -\n^\. (60) 

V 71" UJ 

Let us check now that By{x) Se'Sg for a; ~ at the inflow side of the shock. On one hand, we can expect that 
near the diffusion region 

^ 2e'B'^ = 2Boe (^1 + ^ In ^) (l - f In ^) = 2Bo£ (^1 + ^ In Q . (61) 

On the other hand, for x ~ Id, the global solution tends to 

By ^2Boe {In ^ + 3)] ^ 2Boe ( I + —In^] . (62) 

V ^ Jx=i^ V hJ 

Therefore By{x) — > ^e'Bq, if we take into account the next term in the e expansion for By at the shock. This resolves 
the question concerning the asymptotic behaviour By{x/ld) — *■ 2e' Bd when x/ld — > oo. 

Similarly it can be shown that the global Petschek solution tends to the local one at the distance x ^ Id- This 
implies that all components of V, B are matched automatically near the boundary with the convection zone if one of 
them {Bx in our case) has been adjusted properly. 

Now we can estimate the integral used in equation (|^: 

1 ? (B.(x'.0) <(x',0))^^, ^ 1 t 1 'I ^1 ? 

TT J X' J J n J 

— L — L —l^ la 

The integral over the diffusion region x G {—Id, Id) is estimated as 0{e) since By{x, 0) — By {x, 0) is an odd function 
of X, and the integral converges in the usual sense rather than to be calculated as a principal value. The contribution 
from the intervals {—L,—ld) and {ld,L) are estimated as 0{e^\w{L /Id)) because as it follows from equation ( |6^), the 
difference By{x,0) ~ By{x,0) ~ 0{e^\n{L/ld))- Taking into account the hierarchy of the small parameters (|37D we 
conclude that the whole integral ( |6^ ) is estimated as 0{e). 
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Figure Captions 

Figure 1: Scheme of Pctschek reconnection. 

Figure 2: Numerical results for Petschek-type reconnection with localized resistivity. Left column: structure of 
magnetic field lines (solid lines) and stream lines (dashed), distributions of the Vx,Bx, and convection electric field. 
Right column: distributions of the electric current, Vy,By, and total electric field. 

Figure 3: Numerical results for Sweet Parker reconnection with constant resistivity. Left column: structure of 
magnetic field lines (solid lines) and stream lines (dashed), distributions of the Vx,Bx, and convection electric field. 
Right column: distributions of the electric current, Vy,By, and total electric field. 



12 



This figure "Figurel.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/physics/0 11 1 188vl 



This figure "Figure2.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/physics/0 11 1 188vl 



This figure "Figure3.gif" is available in "gif" format from: 



http://arXiv.0rg/ps/physics/0 11 1 188vl 



